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Abstract 



We report the results of extensive Dynamic Monte Carlo simulations of 
systems of self-assembled Equilibrium Polymers without rings in good sol- 
vent. Confirming recent theoretical predictions, the mean-chain length is 
found to scale as (L) = L* oc cff' exp(5£^) with exponents = 5^ = 

1/(1 + 7) ~ 0.46 and = [1 + (7 - l)/{vd - l)]/2 ^ 0.60, 5^ = 1/2 in the 
dilute and semi-dilute limits respectively. The average size of the micelles, as 
measured by the end-to-end distance and the radius of gyration, follows a very 
similar crossover scaling to that of conventional quenched polymer chains. In 
the semi-dilute regime, the chain size distribution is found to be exponential, 
crossing over to a Schultz-Zimm type distribution in the dilute limit. The 
very large size of our simulations (which involve mean chain lengths up to 
5000, even at high polymer densities) allows also an accurate determination 
of the self-avoiding walk susceptibility exponent 7 = 1.165 ± 0.01. 
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Systems in which polymerization is beheved to take place under condition of chemical 
equilibrium between the polymers and their respective monomers are termed "equilibrium 
polymers" (EP). Examples include liquid sulfur [jl],^ and selenium |p, poly-a-methylstyrene 
0, giant micelles P| and protein filaments 0. In several of these systems, the behaviour is 
strongly affected by the presence of polymeric rings [|^ . For reasons that are not entirely clear 
however, ring-formation effects seem to be negligible in other cases, including that of giant 
micelles [Q. Since the latter are among the most widely studied examples of equilibrium 
polymers, the ring-free case is of some interest. For micelles there has been controversy 
concerning the dependence of the mean chain length on volume fraction, described by the 
growth exponent a. Indirect evidence, based on viscosity measurements interpreted via the 
reptation-reaction model P], suggests in some semidilute systems a value of a ^ 1.2 different 

from the prediction of scaling theory a ~ 0.6 p|-p . 



In this communication we present results for EP without rings, derived by means of a new 
efficient Dynamic Monte Carlo algorithm [jT^ , sampling over a wide range of volume fractions 
and scission energies E. (This is the energy to create a pair of chain ends.) We confirm 
unambiguously the scaling expectations based on the classical behaviour of conventional 
quenched polymers ||9HrT|. The crossover scaling of both mean chain length (L) and mean 



size of the chains are discussed. Our results provide strong support for the form of the 
Molecular Weight Distribution (MWD), c(L), predicted on the ground of renormalizational 
analysis |lTO| , [rT| , and discard an earlier conjecture [0. In passing, we obtain an accurate 



estimation for the susceptibility exponent 7 of self-avoiding walks, which can be directly 
measured from the slope of c(L) in semi-log coordinates. 

At the level of a Flory-Huggins Mean Field Approximation (MFA) the grand potential 
r2 for a system of EP may be written, in units of k^T ^ as 



= ^ c[L)h'^ [in (c(L)6'*) + E + L/i] . (1) 



The first term is the entropy of mixing, the second the scission energy of a bond and the 
third entails the usual Lagrange multiplier /i for the constraint of constant total density of 
the monomers = Y^J^ Lc{L) (in suitable units). For convenience the persistence length 
Ip is set equal to the size of a monomer h. Without loss of generality we have set to zero 
in Eq. (^ the part of the free energy linear in chain length. It is useful to introduce the 
normalized probability distribution p{L) = ((L)/0)c(L), so that J2lP{L) = 1- Minimizing 
with respect to the MWD, subject to the constraint, yields an exponentially decaying MWD, 
with p(L)dL = p{s)ds = exp{~s)ds, where the chemical potential sets the scaling variable 
s = fiL. This distribution requires {L)fi = 1, i.e. the scaling variable s is given, in the 
MFA, by the reduced chain length x = L/ (L) = s. The mean chain length is found to 
be (L) oc 0^/^ exp(£'/2). This result is expected to be a good approximation around the 9- 
temperature of the system (not of interest here) and also in the dense regime where excluded 
volume interactions are largely screened, and the chains behave like (gaussian) random walks 



1^ . The MFA coincides, in the low concentration limit, with the law of mass action. 

One can easily extend the above analysis to dilute and semi-dilute solutions of EP. We 
recall [1^ from ordinary polymers that the correlation length S for a chain of length L in the 



dilute limit is given by the size, R, of the chain: E = R oc L'^. When chains become so long 
that they start to overlap at L ^ L* oc 0"^/'^'^'^"^^ the correlation length of the chain levels 
off and reflects the (chain-length independent) 'blob-size', C,'- thus S = ^ oc L*". Here d is the 



1 



dimension of space and v ~ 0.588 is the swollen chain exponent. The introduction of these 
correlations into the picture is relatively straightforward. The mean-field approach remains 
valid 1^ as long as the basic 'monomer' is replaced by coarse-grained blobs of monomers, of 
size S. Accordingly we may write the grand potential as 

9. = Y, c{L)E'' [in (c(L)S'^) - In {E/bf^' + E + Lfi] . (2) 

L 

In Eq.(^ the second term in the brackets reflects the gain in entropy when a chain breaks 
so that the two new ends can explore a volume Entropy is increased because the 
excluded volume repulsion on scales less than S is reduced by breaking the chain; this effect 
is accounted for by the additional exponent 6 = (7 — l)/z/. Note that in 3d one has 6 ~ 0.3 
with 7 ^ 1.17. In MFA 7 = 1, 6 = 0, and Eq. (H) simplifies to Eq. (|l]), where correlations, 
brought about by the mutual avoidance of the chains, i.e. excluded volume, are ignored. In 
the dilute regime Eq. (0) can be rewritten as = J2l c(-Z^) [In c(L) + (7 — 1) InL + -|- L/i], 
so that the relation to the well known partition function of self-avoiding walks [1^ (with an 
effective coordination number z), Ql oc z^L'^~^, is evident. Hence Eq.(^ covers the entire 
concentration range. 

Minimizing Eq. at fixed 0, yields the MWD p{L)dL = p{s)ds oc H(s) exp{—s)ds. Note 
that for large chain length L we have p{L)dL oc exp(— 7ggx)(ix. Here, we have introduced 
the effective exponent 7gg = s/x = {L)fi = (s). Since in general 7gg is not unity, the 
scaling variable s does not equal the reduced chain length x. In the two limits, far away 
from the crossover line dividing dilute from semidilute behavior, one can readily calculate 
the integrals. In the semi-dilute limit we obtain 7gg = 1, so that in this case s = x as in 
conventional MFA. In the dilute limit, however, 7^^ = 7, so that s = 'jx. Substituting these 
results, the distributions become 

pyx)ax - I _g_3;7-i exp(-7x)rfx ((L) < L*) ^'^^ 
Using again the conservation law one obtains a mean chain length 

(L) = oc 0" exp(5E) (4) 

with exponents ad = 5d = 1/(1 + 7) ~ 0.46 in the dilute and = 1/2(1 + (7 — 1)/ {ud— 1) ^ 
0.6, 5s = 1/2 in the semi-dilute regime. Note that for self-consistency we require (L)^ = {L)^ 
at = (J)*, which imposes that the crossover density obey 0* oc exp{—E/(p) and the crossover 
chain length L* oc exp^E/n). Here the exponents are ip = («<; — ad)/{Ss — Sd) ~ 3.8 and 
K = {ud — l)(f ^ 2.93; similar relations hold between the amplitudes. 

We have confirmed all the major predictions of this scaling analysis by means of an 
efficient new Dynamic Monte Carlo simulation, based on a modified bond fluctuation algo- 
rithm for polymer chains [T^ : in EP systems bonds between monomers break and recombine 



permanently and the chains are only transient objects; the (two saturated or unsaturated) 
bonds of each monomer become the principle objects of the data structure |T^. We were 
able to obtain mean chain lengths of up to (L) ^ 5000 with large system sizes (e.g. we have 
62,500 monomers at = 0.5). These values are large enough for meaningful comparison 
with experimental systems of giant micelles, and represent at least an order of magnitude 
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FIG. 1. Crossover scaling versus reduced density <j)/(l)* for (a) rescaled mean chain length {L)/L* 
confirming scaling Eq. @ and (b) reduced mean-square chain size u (Same symbols as in (a), stars 
for dead chains). 



increase over those studied previously |16[. (They also represent higher mean chain lengths 
at melt densities than any we have so far seen published for conventional quenched chains.) 
Details of the method, and more of our new results, will be presented elsewhere [|12|. Note 



that our simulation is limited by system size rather than equilibration time: we are unable 
to go to higher E (lower temperatures) to increase (L) further because the largest chain 
present would then contain too high a proportion of all the monomers present. This leads 



to finite size corrections, studied in detail in Ref. . 

The growth exponents a for dilute and semidilute regimes are most readily analysed by 
plotting the number of blobs {L)/L* versus the reduced density 0/0*. (This is suggested by 
the scaling prediction Eq. @.) The data are found to collapse onto a single master curve in 
Fig. |I|a with remarkable precision; figure represents one of the main results of this work. The 
two limiting slopes compare very well with the two predicted growth exponents {ad = 0.46, 
as ~ 0.60). In addition, we are able to confirm the other exponents mentioned above and 



to obtain all the corresponding amplitude factors |]T2 

We show in Fig. the crossover scaling plot of the mean chain size as measured from 
the end-to-end distance {Rl) and the radius of gyration (R^^- Defining a reference chain 

size Rq = b{Ly, we plot in Fig. |pD the reduced average chain size, Ue = (Rl) / Rq'^ and 

'2\ / D 2 



Ug = {Rgj/RQ , defined via the end-to-end distance and the radius of gyration respectively, 
as functions of a scaling variable v oc 0/0*. The curves for EP are compared to those of 
conventional quenched monodisperse polymer (so that c(L) = 5{L — N)). The latter data, 
represented as stars in Fig. ^o, was taken from Paul et al |T^. In the semi-dilute limit 
the chains are Gaussian on length scales larger than the blob size ^ implying the scaling 
u oc with exponent (3 = (2z/ — l)/(3z/ — 1) ^ —0.23. It is remarkable that, in their 
scaling behaviour, EP and dead polymers are nearly indistinguishable; the two universal 
scaling functions coincide to numerical accuracy. 

The normalized MWD curves, p{x) plotted versus the reduced chain length x = L/ (L), 
likewise collapse perfectly on single 'master' curve as shown in Fig. ^jfor systems in the dilute 
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limit. With larger values of the reduced chain length x one observes an increasing statistical 
noise since very long chains are more rarely encountered in the simulation. Our results are 
qualitatively consistent with the Schultz distribution Eq. (^d) (bold line) and agree with 
earlier studies of the MWD of EP in two dimensions [T^ where the much larger value of 



l2d = 43/32 strongly enhances the effect. In three dimensions the difference is smaller 
yet detectable. We have also included the high-density prediction so as to emphasize the 
observed deviation. As shown in the insert of Fig. ^a, our MWD in dilute solutions is 
qualitatively consistent with the additional power-law dependence p oc x'^~^ in the limit of 
small X. (Note that maximum of the distribution is at Xmax = (7 ~ l)/7 ~ 0.1 where we 
have chain lengths of only L ^ 4.) 

In the semi-dilute regime (not shown) the normalized MWD shows an equally good 
data collapse as a function of x, confirming the prediction Eq. (^) |]12|. At intermediate 
densities, slightly above the dilute-semi-dilute crossover line, a non-negligible fraction of 
chains are smaller than the blob size, and being thus fully swollen, may exist within the 
network made up from the chains of average size (L) without being seriously perturbed by 
the interchain interaction. The molecular weight distribution crosses therefore smoothly over 
from the dilute limit to the semi-dilute limit; detailed data will be presented elsewhere |]12 



Here we elucidate this crossover by analyzing the effective exponent 7gg found from fitting 
to the expected exp(— 7gga;) tail in the MWD at large x. The obtained values are plotted 
in Fig. 03 versus the number of blobs {L)/L*. It is clearly seen that 7gg — > 7 = 1.165 
and 7gg — 7 = 1 in the dilute and semi-dilute limits respectively. In between we observe 
a smooth crossover for the effective exponent. The error bars are around ±0.01, however. 



a much higher accuracy is in principle feasible with our method |]T2|. Note also that the 
heterogeneity index / = (L^) /{L)"^ = 1 + I/7, (the ratio of weight average and number 
average molecular weights) allows for an alternative determination of 7. However, due to the 
strong contribution of the distribution at small x, which is - as mentioned above - influenced 
by finite chain- length effects, this method systematically overestimates the exponent 

In summary, we have presented new simulation data on the statistics of Equilibrium 
Polymers. By use of an efficient Dynamic Monte Carlo algorithm, we are able to study much 
larger mean chain lengths (combined with high densities) than any previously examined. The 
present work unambiguously confirms the scaling results for EP without rings, based on the 



classical behaviour of conventional quenched polymers P-[Tl[|. In particular our results do not 
support a value of the semidilute growth exponent a larger than the classical one {a ~ 0.6), 
as suggested previously in explanation of the anomalous viscosity of giant micelles in 
nonpolar solvents. We also find no evidence for a possible negative exponent in the power 
law in Eq.(^, as postulated in some treatments of the unusual diffusive behaviour in giant 
micelles ||19[. The extremely large size of our chains, allied to a careful analysis of the size 
distribution in the dilute limit, allows us to extract an accurate estimate of the self- avoiding 
walk susceptiblity exponent 7 = 1.165 ± 0.01. This compares with the best renormalization 
group estimate 7 = 1.1615 ± 0.0011 |2|]. 
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FIG. 2. Molecular Weight Distribution, (a) Dilute limit confirming p{x) oc exp(— 72;) with 
7 ~ 1.165. Insert: MWD for small x for two configurations with (j) = 0.005 (squares) and (f) = 0.01 
(diamonds) both at E = 10 in the dilute regime. The prediction of Eq. (|3|a) is also given for 
comparison (bold line), (b) Effective exponent 7ggf versus the number of blobs per chain {L)/L*. 

burgh (TRACS program). JW is indebted to D.P. Landau for hospitality at the University 
of Georgia, Athens, USA. 
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